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"Metodo di post-elaborazione per ridurre artefatti in immagini digitali 
codificate a blocchi, e dispositivi di post-elaborazione per attuare tale metodo." 

DESCRIZIONE 

La presente invenzione si riferisce alia codifica di immagini digitali. Piu 
precisamente, l'invenzione riguarda un metodo di post-elaborazione per ridurre 
artefatti in immagini digitali codificate a blocchi, e ad un dispositivo di post- 
elaborazione atto ad attuare tale metodo. 

Con la diffusione dei sistemi di comunicazione digitali le immagini digitali 
sono sempre maggiormente impiegate. Cio ha portato alia diffusione di 
telecamere fisse e portatili con capacita di acquisizione ed elaborazione digitali. 

Per meglio sfruttare i mezzi di memorizzazione e la banda trasmissiva 
sono stati sviluppati standard di compressione delle immagini, quali il JPEG 
per le immagini fisse, e MPEG-1 ed MPEG-2 per sequenze di immagini 
televisive. 

I summenzionati standard di compressione prevedono la codifica a blocchi 
basata sulla trasformata coseno discreta (DCT). Un'immagine digitale viene 
suddivisa in blocchi di pixel, e ciascun blocco viene codificato 
indipendentemente dagli altri. Vengono calcolati coefficienti DCT per i pixel di 
ciascun blocco, e ai coefficienti DCT viene applicata una matrice di 
quantizzazione per ridurre l'informazione da memorizzare o trasmettere. 
Quando l'immagine deve essere visualizzata, essa deve essere preventivamente 
decodificata. 

A causa del processo di quantizzazione, questi metodi di compressione 
delle immagini danno origine a perdite di informazione neH'immagine 
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decodificata rispetto all'immagine originale. L'immagine decodificata puo 
pertanto presentare un notevole degrado, consistente soprattutto in due tipi di 
artefatti noti nella tecnica con i nomi di "rumore a griglia" ("grid noise") e 
"rumore a gradini" ("staircase noise"). 

Per ridurre il degrado deirimmagine, sono stati proposti metodi di post- 
elaborazione deirimmagine decodificata che consentono di attenuare il rumore 
a griglia ed il rumore a gradini. 

In vista dello stato della tecnica descritto, e uno scopo della presente 
invenzione quello di fomire un nuovo metodo di post-elaborazione per ridurre 
gli artefatti in immagini digitali codificate a blocchi. 

In accordo con la presente invenzione, tale scopo viene raggiunto mediante 
un metodo di post-elaborazione per ridurre gli artefatti in immagini digitali 
codificate a blocchi, caratterizzato dal fatto di prevedere: 

a) la suddivisione di un immagine in ingresso in una pluralita di blocchi di 
immagine; 

b) per ciascun blocco di immagine, la stima di caratteristiche globali di 
detto blocco di immagine che forniscono informazioni sul contenuto medio di 
bordi di immagine lungo le direzioni orizzontale e verticale di detto blocco di 
immagine; 

c) per ciascun pixel di un blocco di immagine in esame, la stima di 
caratteristiche locali relative a detto pixel che forniscono informazioni sul 
contenuto di bordi di immagine lungo la direzione orizzontale e verticale di 
un'area di immagine neH'intorno di detto pixel; 

d) la modifica del valore di detto pixel in funzione sia di dette 
caratteristiche globali del blocco di immagine cui appartiene detto pixel, sia 
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delle caratteristiche locali dell'area di immagine nell'intomo di detto pixel. 

Sempre in accordo con l'invenzione, viene fornito un dispositivo di post- 
elaborazione per attuare tale metodo, caratterizzato dal fatto di comprendere: 

- primi mezzi alimentati da un' immagine di ingresso per stimare 
caratteristiche globali di un blocco di immagine in esame, dette caratteristiche 
globali fomendo informazioni su un contenuto medio di bordi di immagine 
lungo la direzione orizzontale e verticale di detto blocco di immagine; 

- secondi mezzi alimentati da detta immagine di ingresso per stimare 
caratteristiche locali per ciascun pixel del blocco di immagine in esame, detta 
caratteristiche locali fornendo informazioni sul contenuto di bordi di immagine 
lungo le direzioni orizzontale e verticale di un'area di immagine neH'intorno di 
detto pixel; 

- terzi mezzi alimentati da dette caratteristiche globali e da dette 
caratteristiche locali per modificare il valore di detto pixel in funzione sia di 
dette caratteristiche globali che di dette caratteristiche locali. 

Le caratteristiche ed i vantaggi della presente invenzione saranno resi 
maggiormente evidenti dalla seguente descrizione dettagliata di una sua forma 
di realizzazione, descritta a titolo di esempio non limitativo negli uniti disegni, 
nei quali: 

la figura 1 e un diagramma a blocchi schematico che illustra il principio di 
funzionamento di un metodo secondo la presente invenzione; 

la figura 2 mostra un' immagine digitale suddivisa in blocchi di immagine; 

la figura 3 mostra in dettaglio uno dei blocchi di immagine in cui e 
suddivisa rimmagine di figura 2; 

la figura 4 mostra una matrice di pixel del blocco di immagine di figura 3; 
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la figura 5 mostra un sottoblocco di immagine del blocco di immagine di 
figura 3 utilizzato valutare le caratteristiche globali del blocco di immagine; 

la figura 6 mostra una finestra di elaborazione orizzontale utilizzata per 
valutare caratteristiche locali nella direzione orizzontale per un generico pixel 
del blocco di immagine; 

la figura 7 mostra una finestra di elaborazione verticale utilizzata per 
valutare caratteristiche locali nella direzione verticale per detto generico pixel; 

le figure 8 e 9 mostrano due funzioni di appartenenza utilizzate per 
eseguire un calcolo fuzzy; 

la figura 10 e un diagramma a blocchi di un dispositivo secondo la 
presente invenzione; 

la figura 1 1 mostra la struttura di due blocchi del dispositivo di figura 10; e 

la figura 12 e Un diagramma a blocchi di altri due blocchi del dispositivo di 
figura 10. 

Con riferimento alia figura 1, vi e mostrato un diagramma a blocchi che 
illustra il principio di funzionamento del metodo di post-elaborazione secondo 
la presente invenzione. Un'immagine digitale di ingresso compressa e 
decodificata In viene alimentata ad un blocco di estrazione delle caratteristiche 
FE. II blocco FE analizza rimmagine per valutare le sue caratteristiche globali 
e locali. Le caratteristiche globali e locali, rispettivamente GF e LF, 
deirimmagine In valutate dal blocco FE sono alimentate in ingresso ad un 
blocco di elaborazione Fuzzy FUZZY che, in accordo a regole fuzzy, determina 
parametri FA atti a determinare il tipo di filtraggio da eseguire, in accordo alle 
caratteristiche globali e locali GF e LF deirimmagine. I parametri FA calcolati 
dal blocco FUZZY sono alimentati ad un blocco di composizione del filtro FC 
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che, in accordo a detti parametri FA, determina il tipo di filtraggio da eseguire 
aH'interno di un insieme di filtri predefmiti (blocco FS). II blocco FC 
determina parametri di filtraggio FP che vengono alimentati ad un blocco di 
elaborazione PROC, alimentato direttamente daH'immagine di ingresso In, il 
quale effettua il filtraggio deH'immagine di ingresso In in accordo ai parametri 
di filtraggio FP, per fomire una immagine post-elaborata di uscita Out. 

Si evince che il tipo di filtraggio da eseguire sull' immagine decodificata di 
ingresso In viene scelto dopo una stima delle caratteristiche globali e locali 
deirimmagine decodificata di ingresso. Per le area dell'immagine vicine a 
rumore a griglia e nei pressi di un bordo viene eseguito un filtraggio passa- 
basso, in modo da ridurre sia il rumore a gradini che il rumore a griglia. Per 
quelle aree contenenti dettagli di piccole dimensioni (bordi deH'immagine) non 
viene effettuato alcun filtraggio, Pertanto, il metodo secondo la presente 
invenzione prevede di eseguire un filtraggio non lineare adattativo sui pixel 
deirimmagine decodificata. 

11 principio di funzionamento sopra delineato verra ora descritto in 
dettaglio. 

Come mostrato in figura 2, 1' immagine di ingresso In e suddivisa in 
blocchi di immagine EB, ciascuno contenente un ugual numero di pixel. Una 
dimensione tipica dei blocchi e 8*8 pixel (figura 3), ma cio non deve tuttavia 
essere inteso come una limitazione, poiche sono possibili altre dimensioni per i 
blocchi. 

I blocchi di immagine IB deirimmagine di ingresso In vengono scanditi 
linea per linea a partire dal blocco in alto a sinistra fino a quello in basso a 
destra. II blocco di estrazione delle caratteristiche FE di figura 1 determina per 
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ciascun blocco IB le caratteristiche globali e locali GF, LF. 

Le caratteristiche globali del blocco di immagine IB in esame vengono 

determinate applicando operatori di Sobel orizzontali e verticali: 

"hll hl2 hl3~ 
operatore di Sobel orizzontale (Hsob): h21 h22 h33 

h31 h32 h33 



operatore di Sobel verticale (Vsob): 



vll vl2 vl3 
v21 v22 v23 
v31 v32 v33 



a ciascun pixel appartenente ad un sotto-blocco di immagine interno al 

blocco di immagine IB. Ad esempio, i seguenti operatori di Sobel: 
-1 0 1 



Hsob: 



2 0 2 
-1 0 1 



Vsob: 



-1 
0 
1 



-2 -1 
0 0 
2 1 



vengono applicati ad un sotto-blocco di immagine ISB di 6*6 pixel (area 
grigia in figura 5). Come mostrato in figura 4, per ciascun pixel Pi del sotto- 
blocco di immagine ISB viene presa in considerazione una matrice M di 3*3 
pixel vicini centrata nel pixel Pi, ed i valori dei pixel di detta matrice vengono 
moltiplicati per i coefficienti degli operatori di Sobel orizzontali e verticali, per 
ottenere: 

Hsob = (P3+P8+2*P5)-(P1+P6+2*P4), 
Vsob=(P6+P8+2*P7)-(P 1 +P3+2*P2), 

ove P1-P4 e P5-P8 sono i valori dei pixel (livelli di grigio). 



6 



Come noto, gli operatori di Sobel orizzontali e verticali effettuano un 
filtraggio in grado di determinare bordi lungo la direzione orizzontale e 
verticale, rispettivamente. 

I valori di uscita degli operatori di Sobel orizzontali calcolati per i pixel del 
sotto-blocco di immagine ISB vengono accumulati per ottenere un valore 
accumulato Acc(Hsob), ed i valori di uscita degli operatori di Sobel verticali 
calcolati per i pixel del sotto-blocco di immagine ISB sono accumulati per 
ottenere un valore accumulato Acc(Vsob). Acc(Hsob) fomisce il contenuto di 
alte frequenze nella direzione orizzontale (bordi verticali) del blocco di 
immagine EB. Acc(Vsob) fomisce il contenuto di alte frequenze nella direzione 
verticale (bordi orizzontali) del blocco di immagine IB. Pertanto, Acc(Hsob) e 
Acc(Vsob) forniscono rispettivamente il grado di presenza di bordi del blocco 
di immagine in esame nelle direzioni verticale ed orizzontale. Si deve notare 
che per valutare le caratteristiche globali GF del blocco di immagine IB in 
esame solo i pixel appartenenti a questo blocco vengono presi in 
considerazione (applicando gli operatori di Sobel 3*3 al sotto-blocco di 
immagine ISB, non e necessario prendere in considerazione i pixel appartenenti 
ai blocchi di immagine vicini). 

Le caratteristiche globali GF del blocco di immagine in esame possono 
essere costituite dai valori accumulati Acc(Hsob) e Acc(Vsob). In alternativa, 
le caratteristiche globali GF del blocco di immagine possono essere costituite 
da una media dei valori accumulati Acc(Hsob) e Acc(Vsob), per ottenere il 
numero medio di bordi nelle direzioni orizzontale e verticale del blocco di 
immagine in esame. 

Le caratteristiche locali del blocco di immagine EB vengono stimate 
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secondo il metodo seguente. Tutti i pixel del blocco di immagine IB in esame 
vengono scanditi linea per linea partendo dal pixel in alto a sinistra fino a 
quello in basso a destra. A ciascun pixel e applicata una finestra di elaborazione 
orizzontale: viene preso in considerazione un numero prestabilito di pixel 
appartenenti alia stessa linea di immagine del pixel in considerazione e che 
precedono e seguono tale pixel. Una finestra di elaborazione orizzontale HPW 
adatta alio scopo, ossia una finestra di elaborazione orizzontale 1*5, e mostrata 
in figura 6: per un dato pixel, si considerano i due pixel che precedono Pa, Pb 
ed i due pixel che seguono Pc, Pd appartenenti alia stessa linea. In figura 6 
viene mostrata ad esempio la finestra di elaborazione HPW associata al primo 
pixel Px del blocco di immagine. Si deve notare che non solo vengono 
considerati i pixel del blocco di immagine EB in esame, ma anche i pixel 
appartenenti a blocchi di immagine vicini; cio e ad esempio il caso del primo, 
del secondo, del settimo e dell'ottavo pixel di ciascuna linea di pixel del blocco 
di immagine IB in esame. 

L'operatore di Sobel orizzontale summenzionato viene applicato a ciascun 
pixel Pa, Pb, Px, Pc, Pd della finestra di elaborazione orizzontale HPW, per 
ottenere cinque valori di uscita HS1-HS5. I valori HS1-HS5 fomiscono le 
caratteristiche locali nella direzione orizzontale per il pixel in esame Px, ossia il 
contenuto di alte frequenze nella direzione orizzontale della regione 
deirimmagine neH'intorno del pixel in esame. 

Analogamente, a ciascun pixel del blocco di immagine IB viene applicata 
una finestra di elaborazione verticale. La finestra di elaborazione verticale e 
costituita dal pixel in esame Px e da un numero prestabilito di pixel 
appartenenti alia stessa colonna e che precedono e seguono il pixel in esame; 
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ad esempio, come mostrato in figura 7, la finestra di elaborazione verticale 
VPW puo avere dimensioni identiche alia finestra di elaborazione orizzontale 
HPW (5*1), e contenere quindi due pixel Pe, Pf che precedono il pixel Px e due 
pixel Pg, Ph che seguono il pixel Px nella direzione verticale. 

L'operatore di Sobel verticale Vsob precedentemente descritto viene 
quindi applicato a ciascun pixel della finestra di elaborazione verticale VPW 
per ottenere cinque valori VS1-VS5. I valori VS1-VS5 formano le 
caratteristiche locali nella direzione verticale per il pixel in esame, ossia il 
contenuto di alte frequenze nella direzione verticale di una regione di 
immagine neH'intomo del pixel in esame. 

Le caratteristiche globali GF del blocco di immagine IB in esame (ossia i 
due valori accumulati Acc(Hsob) e Acc(Vsob) o, in alternativa, il valor medio 
di Acc(Hsob) e Acc(Vsob)) e le caratteristiche locali LF per il pixel in esame 
all'interno di detto blocco di immagine (i dieci valori HS1-HS5 e VS1-VS5) 
vengono quindi alimentati al blocco di calcolo fuzzy FUZZY. Questo blocco 
valuta il gradi di appartenenza di un generico valore HSi e VSi (i=1..5) a due 
insiemi fuzzy "Small" ("Piccolo") e "Big" ("Grande"); questi gradi di 
appartenenza possono essere valutati applicando a HSi, VSi le funzioni di 
appartenenza raffigurate nelle figure 8 e 9. In queste figure, Thl e Th2 sono 
valori dipendenti dalla caratteristiche globali GF del blocco di immagine in 
esame, ossia dai valori accumulati Acc(Hsob) e Acc(Vsob) o dal valor medio 
dei valori accumulati; nel primo caso, Thl e Th2 sono diversi per i valori HSi e 
VSi; nel secondo caso, Thl e Th2 sono gli stessi per i valori HSi e VSi. 

Vengono poi valutate regole fuzzy aventi come antecedenti io gradi di 
appartenenza dei valori di uscita HSi e VSi ai due insiemi fuzzy "Piccolo" e 
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"Grande". Cio significa che devono essere valutate 32 regole sia per la 
direzione orizzontale che per quella verticale. Tuttavia, tutte le regole fuzzy 
aventi la stessa conseguenza vengono sintetizzate in una regola soltanto 
mediante un operatore Else. In questo modo, la complessita del sistema viene 
ridotta, e deve essere calcolato un totale di nove regole per ciascuna direzione. 

Le seguenti regole fuzzy vengono applicate ai cinque valori HS1-HS5 
associati alia direzione orizzontale: 

1. Se HS1 e Piccolo e HS2 e Piccolo e HS3 e Piccolo e HS4 e Piccolo e 
HS5 e Piccolo, allora ale Grande; 

2. Se HS1 e Piccolo e HS2 e Piccolo e HS3 e Piccolo e HS4 e Piccolo e 
HS5 e Grande, allora cc2 e Grande; 

3. Se HS1 e Piccolo e HS2 e Piccolo e HS3 e Piccolo e HS4 e Grande e 
HS5 e Piccolo, allora cc3 e Grande; 

4. Se HS1 e Piccolo e HS2 e Piccolo e HS3 e Piccolo e HS4 e Grande e 
HS5 e Grande, allora a4 e Grande; 

5. Se HS1 e Piccolo e HS2 e Grande e HS3 e Piccolo e HS4 e Piccolo e 
HS5 e Piccolo, allora a5 e Grande; 

6. Se HS1 e Grande e HS2 e Piccolo e HS3 e Piccolo e HS4 e Piccolo e 
HS5 e Piccolo, allora ct6 e Grande; 

7. Se HS1 e Grande e HS2 e Piccolo e HS3 e Piccolo e HS4 e Piccolo e 
HS5 e Grande, allora ct7 e Grande; 

8. Se HS1 e Grande e HS2 e Grande e HS3 e Piccolo e HS4 e Piccolo e 
HS5 e Piccolo, allora a8 e Grande. 

II livello di attivazione di ciascuna regola dipende dal grado di 
appartenenza dello schema dei valori di uscita HSi dell'operatore di Sobel 
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orizzontale applicato ai cinque pixel della finestra di elaborazione orizzontale 
HPW; il grado di appartenenza dipende a sua volta dalle caratteristiche globali 
GF del blocco di immagine al quale appartiene il pixel in esame. II livello di 
attivazione della regola else (nona regola) e calcolato come aelse = (1 - aave), 
ove aave e la media dei gradi di attivazione delle regole fuzzy da 1 a 8. al-a8 
e aelse, ed un analogo insieme di nove gradi di attivazione delle regole fuzzy 
applicate ai valori VS1-VS5) costituiscono l'uscita FA del blocco di calcolo 
fuzzy FUZZY di figura 1 . 

Ciascuna delle regole sopra elencate e associata ad un rispettivo insieme di 
parametri di filtro predefiniti, che sono memorizzati come una tabella di 
riferimento nel blocco FS di figura 1. Insiemi di parametri di filtro predefiniti 
sono ad esempio: 

Regola 1: (cl 1=1.0, cl2=1.0, cl3=1.0, cl4=1.0, cl5=1.0) se il pixel in esame e 

all'estemo del sotto-blocco di immagine ISB, e 

(cl 1=0.0, cl2=1.0, cl3=1.0, cl4=1.0, cl5=0.0) se il pixel in esame e 

all'interno del sotto-blocco di immagine ISB; 
Regola 2: (c21=0.5, c22=1.0, c23=1.0, c24=1.0, c25=0.0); 
Regola 3: (c31=0.5, c32=1.0, c33=1.0, c34=0.0, c35=0.0); 
Regola 4: (c41=0.5, c42=1.0, c43=1.0, c44=0.0, c45=0.0); 
Regola 5: (c5 1=0.0, c52=0.0, c53=1.0, c54=1.0, c55=0.5); 
Regola 6: (c61=0.0, c62=1.0, c63=1.0, c64=1.0, c65=0.5); 
Regola 7: (c71=0.0, c72=1.0, c73=1.0, c74=1.0, c75=0.0); 
Regola 8: (c81=0.0, c82=0.0, c83=1.0, c84=1.0, c85=0.5); 
Regola else: (c9 1=0.0, c92=0.0, c93=1.0, c94=0.0, c95=0.0). 

I parametri FP del filtro da applicare ai cinque pixel della finestra di 
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elaborazione orizzontale HPW sono calcolati come media pesata dei nove filtri 
sopra descritti, con fattori di peso costituiti dai gradi di attivazione al-a8 e 
aelse delle rispettive regole fuzzy. 

Essendo ai il grado di attivazione della regola i-esima (i=1..9), la nona 
regola fuzzy essendo la regola fuzzy else (a9 = aelse), e cij i coefficient! dell'i- 
esimo filtro (i=1..9, j=1..5), il fattore peso applicato al filtro i-esimo, associato 
alia i-esima regola fuzzy e: 

Fi = ai • cij 

ed i coefficienti Hi del filtro orizzontale finale da applicare ai pixel della 

finestra di elaborazione orizzontale HPW sono dati da: 
9 

cti • cij 

Hj = ^i 

N 

ove N e un fattore di normalizzazione. 

II valore filtrato orizzontalmente Px del pixel Px in esame (al centro delle 
finestra di elaborazione orizzontale) viene poi calcolato come media pesata dei 
valori dei pixel Pa, Pb, Px, Pc e Pd appartenenti alia finestra di elaborazione 
orizzontale HPW, con fattori peso costituiti dai coefficienti Hj: 

Px = Hl*Pa+ H2*Pb+ H3*Px+ H4*Pc+ H5*Pd. 

Analoghi calcoli vengono eseguiti per la direzione verticale, a partire dai 
valori di uscita VS1-VS5 degli operatori di Sobel verticali applicati ai pixel Pe, 
Pf, Px, Pg e Ph della finestra di elaborazione verticale VPW. I coefficienti Vj 
(J=1..5) del filtro per la direzione verticale sono calcolati in modo del tutto 
analogo a quello usato per determinare i coefficienti Hj: 
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9 

Vj = -i=! 

N 



ove pi (i=1..9) sono i gradi di attivazione di nove regole fuzzy per la 
direzione verticale (analoghi quelli sopra elencati per la direzione orizzontale) e 
cij (j=1..5) sono ora i parametri di filtro predefiniti associati alia i-esima regola 
fuzzy per la direzione verticale. I coefficienti Vj vengono poi applicati ai pixel 
nella finestra di elaborazione verticale VPW per calcolare una media pesata 
degli stessi. II valore filtrato del pixel Px in esame, filtrato sia nella direzione 
orizzontale che in quella verticale, viene fornito all'uscita del blocco PROC. 

II valore del pixel Px in esame da moltiplicare per il coefficiente di filtro 
verticale V3 puo essere il valore Px ottenuto dopo aver applicato ai pixel nella 
finestra di elaborazione orizzontale HPW il filtro orizzontale Hj (j=1..5): 
Out = Vl*Pe + V2*Pf + V3*Px + V4*Pg + V5*Ph. 

In alternativa, e possibile valutare dapprima il valore filtrato verticalmente 
Px del pixel in esame: 

Ex = Vl*Pe + V2*Pf + V3*Px + V4*Pg + V5*Ph, 

e poi effettuare il filtraggio nella direzione orizzontale applicando a questo 
valore il rispettivo coefficiente H3 del filtro orizzontale Hj: 

Out = Hl*Pa + H2*Pb + H3*Ex + H4*Pc + H5*Pd. 

La sequenza non ha importanza, la cosa importante da sottolineare essendo 
che al termine del processo il valore del pixel in esame e il risultato di un 
filtraggio sia orizzontale che verticale. 

La figura 10 e uno schema a blocchi di un dispositivo atto ad implementare 
il metodo precedentemente descritto. II dispositivo comprende due blocchi 
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principali: un primo blocco 1 valuta le caratteristiche globali GF dei blocchi di 
immagine IB in cui rimmagine da post-elaborare e suddivisa, ed un secondo 
blocco 2 valuta le caratteristiche locali LF dei pixel dell' immagine ed effettua il 
filtraggio in accordo alle caratteristiche sia locali che globali. 

Si supponga che rimmagine da post-elaborare venga scandita linea per 
linea in ordine sequenziale. II segnale In e un flusso di pixel deH'immagine di 
ingresso scandita linea per linea. II blocco 1 e alimentato dal segnale In; il 
segnale In alimenta anche una cascata di due memorie di linea LM1 e LM2 le 
cui uscite alimentano il blocco 1 . 

All'intemo del blocco 1, il segnale In e le uscite delle memorie di linea 
LM1 e LM2 alimentano un blocco 3 di ritardi di pixel atto ad implementare la 
finestra di pixel 3*3 che viene utilizzata per il calcolo degli operatori di Sobel 
orizzontale e verticale per i pixel del sotto-blocco di immagine 6*6 ISB 
all'mtemo del blocco di immagine IB. II blocco 3 alimenta un blocco 4 che 
calcola le uscite Hsob e Vsob degli operatori di Sobel orizzontale e verticale 
per quei pixel della linea di immagine corrente che appartengono ai sotto- 
blocchi di immagine 6*6 ISB di ciascun blocco di immagine IB. Le uscite 
Hsob e Vsob del blocco 4 vengono alimentate ad un blocco accumulatore 5 nel 
quale esse vengono accumulate. Dopo otto linee di immagine, ossia dopo che e 
stata scandita una linea di blocchi di immagine IB, i valori accumulati 
Acc(Hsob) e Acc(Vsob) (o in alternativa la loro media) per ciascun blocco di 
immagine IB vengono memorizzati in un blocco di memoria 6. 

L'uscita della memoria di linea LM2 alimenta una cascata di altre otto 
memorie di linea LM3-LM10. II blocco 2 e alimentato in parallelo dalle uscite 
delle memorie di linea LM4-LM10. In questo modo, la valutazione delle 
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caratteristiche locali ed il calcolo dei parametri di filtro inizia dopo che il 
blocco \1 ha stimato le caratteristiche globali GF per una linea di blocchi di 
immagine IB. 

AH'intemo del blocco 2, un blocco 7 di ritardi di pixel e alimentato dalle 
uscite di memorie di linea LM4-LM10; mediante le memorie di linea LM4- 
LM10 ed il blocco di ritardo di pixel 7 e possibile implementare la finestra di 
elaborazione verticale 5*1 VPW. Le uscite L4-L10 del blocco di ritardi di pixel 
7 alimentano un blocco 8 che applica l'operatore di Sobel verticale a ciascun 
pixel aH'interno della finestra di elaborazione verticale VPW. Per evitare 
l'impiego di ulteriori memorie di linea, viene preferito un approccio parallelo 
che prevede di calcolare cinque operatori di Sobel verticali in parallelo; le 
uscite dei cinque operatori di Sobel verticali VS1-VS5 vengono alimentate ad 
un filtro fuzzy verticale 9, che e anche alimentato dalle uscite L6-L10 del 
blocco di ritardo di pixel 7 e dall'uscita MOUT della memoria MEM del 
blocco 1 . MOUT fornisce le caratteristiche globali GF del blocco di immagine 
IB correntemente elaborato dal blocco 2, ossia il valore accumulato Acc(Vsob) 
o, in alternativa, la media di Acc(Vsob) e Acc(Hsob). II blocco di filtro fuzzy 
verticale 9 valuta il grado di appartenenza dei valori VS1-VS5 agli insiemi 
fuzzy "Small" ("Piccolo") e "Big" ("Grande") tenendo conto delle 
caratteristiche globali fornite da MOUT, valuta i livelli di attivazione delle 
nove regole fuzzy per la direzione verticale, calcola i coefficienti Vj (j=1..5) del 
filtro verticale ed applica i coefficienti del filtro verticale Vj ai cinque pixel Pe, 
Pf, Px, Pg, Ph nella finestra di elaborazione verticale VPW, per calcolare il 
valore filtrato verticalmente Ex. del pixel al centro della finestra di elaborazione 
verticale. L'uscita del blocco di filtro fuzzy verticale 9 forma il valore filtrato 
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verticalmente Ed del pixel Pd nella fmestra di elaborazione orizzontale HPW 
mostrata in figura 6 ed alimenta direttamente un blocco di filtro fuzzy 
orizzontale 10; l'uscita Pd del blocco 9 alimenta anche una cascata di quattro 
ritardi di pixel D le cui uscite formano rispettivamente i valori filtrati 
verticalmente Ec., Ex., Eh, E§ dei pixel Pc, Px, Pb, Pa nella fmestra di 
elaborazione orizzontale HPW ed alimentano il blocco di filtro fuzzy 
orizzontale 10. 

In parallelo ai blocchi 8 e 9, le uscite L7-L9 del blocco di ritardo di pixel 7 
alimentano un blocco 1 1 che applica ai pixel della fmestra di elaborazione 
orizzontale HPW gli operatori di Sobel orizzontali. Diversamente dagli 
operatori di Sobel verticali, viene calcolato un solo operatore di Sobel 
orizzontale alia volta; un blocco di ritardo di compensazione 12 introduce un 
ritardo per compensare il ritardo di elaborazione del blocco di filtro fuzzy 
verticale 9. L'uscita del blocco 12, che costituisce l'uscita HS5 dell'operatore 
di Sobel orizzontale applicata al pixel Pd della fmestra di elaborazione 
orizzontale di figura 6, alimenta il blocco di filtro fuzzy orizzontale 10 ed una 
cascata di quattro ritardi di pixel D, le uscite dei quali costituiscono i valori 
HS4, HS3, HS2 e HS1 ed alimentano il blocco di filtro fuzzy orizzontale 10. II 
blocco di filtro fuzzy orizzontale 10, che e inoltre alimentato dalFuscita MOUT 
del blocco di memoria MEM nel blocco 1 che fornisce il valore Acc(Hsob) (o 
in altemativa la media dei valori Acc(Hsob) e Acc(Vsob)), valuta il grado di 
appartenenza dei valori HS1-HS5 agli insiemi fuzzy "Piccolo" e "Grande" in 
funzione del valore delle caratteristiche globali GF fomite da MOUT, valuta i 
livelli di attivazione delle nove regole fuzzy precedentemente descritte per il 
filtraggio nella direzione orizzontale, calcola i coefficienti Hj del filtro 
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orizzontale ed applica i parametri Hj ai valori filtrati verticalmente Pa, Pb, Px, 
Ee, Ed dei pixel Pa, Pb, Px, Pc, Pd della fmestra di elaborazione orizzontale 
HPW per ottenere il valore filtrato orizzontalmente e verticalmente Out del 
pixel Px in esame. 

Un circuito di controllo CTRL governa il funzionamento dei blocchi 1, 2 e 
delle memorie di linea LM1-LM10. 

La figura 1 1 mostra la struttura dei blocchi operatori di Sobel orizzontale e 
verticale 8 e 11 di figura 10. Essi sono costituiti in modo noto mediante 
sommatori. 

La figura 12 mostra la struttura sia del blocco di filtro fuzzy verticale 9 che 
del blocco di filtro fuzzy orizzontale 10. X1-X5 sono rispettivamente le uscite 
dell'operatore di Sobel verticale e orizzontale VS1-VS5 e HS1-HS5. X1-X5 
vengono alimentate ad un blocco di valutazione delle regole 13 che valuta i 
gradi di attivazione P1-P9 delle nove regole fuzzy per la direzione verticale o, 
rispettivamente, i gradi di attivazione ccl-a9 delle nove regole fuzzy per la 
direzione orizzontale. I gradi di attivazione valutati dal blocco 13 sono 
alimentati ad una tabella di riferimento contenente rispettivi parametri di filtro 
predefmiti F1-F9 (costituente il blocco FS in figura 1), e le uscite della tabella 
di riferimento, ossia i parametri di filtro predefmiti cij moltiplicati per il grado 
di attivazione della rispettiva regola fuzzy, vengono alimentati ad un blocco di 
composizione del filtro 14 che calcol.a i coefficienti VI -V5 o, rispettivamente, 
H1-H5, del filtro verticale o rispettivamente del filtro orizzontale. Detti 
coefficienti vengono poi alimentati ad un blocco di elaborazione 1 5 che e anche 
alimentato dai valori dei pixel PXS (L6-L10 o, rispettivamente, Pa, P_b_, Ex, Ec, 
Pd di figura 10); il blocco 15 applica i coefficienti di filtraggio ai valori dei 
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pixel per ottenere il valore filtrato del pixel in esame Px. di figura 10); il blocco 
15 applica i coefficienti di filtraggio ai valori dei pixel per ottenere il valore 
filtrato del pixel in esame Px. 
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RIVENDICAZIONT 

1 . Metodo di post-elaborazione per ridurre gli artefatti in immagini digitali 
codificate a blocchi, caratterizzato dal fatto di prevedere: 

a) la suddivisione di un immagine in ingresso (In) in una pluralita di 
blocchi di immagine (IB); 

b) per ciascun blocco di immagine (IB), la stima di caratteristiche globali 
(GF) di detto blocco di immagine che fomiscono informazioni sul contenuto 
medio di bordi di immagine lungo le direzioni orizzontale e verticale di detto 
blocco di immagine (IB); 

c) per ciascun pixel (Px) di un blocco di immagine (IB) in esame, la stima 
di caratteristiche locali (LF) relative a detto pixel (Px) che fomiscono 
informazioni sul contenuto di bordi di immagine lungo la direzione orizzontale 
e verticale di un'area di immagine nell'intorno di detto pixel (Px); 

d) la modifica del valore di detto pixel (Px) in funzione sia di dette 
caratteristiche globali (GF) del blocco di immagine (IB) cui appartiene detto 
pixel (Px), sia delle caratteristiche locali (LF) dell'area di immagine 
nell'intorno di detto pixel (Px). 

2. Metodo secondo la rivendicazione 1, caratterizzato dal fatto che la fase 
d) prevede di: 

dl) defmire un insieme di caratteristiche locali predefinite; 

d2) defmire gradi di coincidenza (al-a9,(3l-(39) di dette caratteristiche 
locali (LF) dell'area di immagine neH'intomo di detto pixel (Px) con ciascuna 
delle caratteristiche locali predefinite di detto insieme, detti gradi di 
coincidenza (ccl-a9,pl-P9) essendo dipendenti da dette caratteristiche globali 
(GF) del blocco di immagine (IB) cui appartiene detto pixel (Px); 
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d3) rendere il valore di detto pixel (Px) pari ad una media pesata ( Px . Px ) 
del valore di detto pixel (Px) e dei valori di pixel vicini (Pa-Pd,Pe-Pf), con 
fattori peso (H1-H5,V1-V5) funzione di detti gradi di coincidenza di dette 
caratteristiche locali (LF) con ciascuna di dette caratteristiche locali predefmite. 

3. Metodo secondo la rivendicazione 2, caratterizzato dal fatto che detto 
determinare i gradi di coincidenza (al-cc9,pi-P9) nella fase d2) prevede di 
effettuare un calcolo fuzzy. 

4. Metodo secondo la rivendicazione 3, caratterizzato dal fatto che 
ciascuna di dette caratteristiche locali predefmite e associata ad un rispettivo 
gruppo di fattori peso predefmiti (cij), e che ciascuno di detti fattori peso (Hl- 
H5,V1-V5) e calcolato come media pesata di corrispondenti fattori peso 
predefmiti (cij) di detti gruppi con coefficienti peso (al-a9,pi-{39) funzione di 
detti gradi di coincidenza di dette caratteristiche locali (LF) con ciascuna di 
dette caratteristiche locali predefmite. 

5. Metodo secondo la rivendicazione 4, caratterizzato dal fatto che detta 
stima delle caratteristiche globali (GF) del blocco di immagine (IB) in esame 
prevede di applicare operatori di Sobel orizzontali e verticali (Hsob,Vsob) a 
pixel appartenenti ad un sotto-blocco di immagine (ISB) intemo a detto blocco 
di immagine (IB) in esame. 

6. Metodo secondo la rivendicazione 5, caratterizzato dal fatto che detta 
stima delle caratteristiche globali (GF) del blocco di immagine (IB) in esame 
prevede di sornmare uscite degli operatori di Sobel orizzontali (Hsob) applicati 
a ciascun pixel del sotto-blocco di immagine (ISB) per ottenere una uscita 
accumulata degli operatori di Sobel orizzontali (Acc(Hsob)), e di sommare 
uscite degli operatori di Sobel verticali (Vsob) applicati a ciascun pixel di detto 
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sotto-blocco di immagine (ISB) per ottenere una uscita accumulata degli 
operatori di Sobel verticali (Acc(Vsob)). 

7. Metodo secondo la rivendicazione 6, caratterizzato dal fatto che dette 
caratteristiche globali (GF) del blocco di immagine (IB) in esame sono 
costituite da dette uscite accumulate degli operatori di Sobel orizzontali e 
verticali (Acc(Hsob),Acc(Vsob)). 

8. Metodo secondo la rivendicazione 6, caratterizzato dal fatto che dette 
caratteristiche globali (GF) del blocco di immagine (IB) in esame sono 
costituite da una media di dette uscite accumulate degli operatori di Sobel 
orizzontali e verticali (Acc(Hsob),Acc(Vsob)). 

9. Metodo secondo la rivendicazione 7 o 8, caratterizzato dal fatto che 
detta stima delle caratteristiche locali (LF) per un pixel (Px) del blocco di 
immagine (IB) in esame prevede di: 

cl) considerare una finestra di elaborazione orizzontale (HPW) contenente 
il pixel (Px) in esame e pixel vicini (Pa,Pb,Pc,Pd) appartenenti ad una stessa 
linea di immagine del pixel (Px) e precedenti e seguenti il pixel (Px); 

c2) applicare detto operatore di Sobel orizzontale (Hsob) a ciascun pixel 
(Pa-Pd,Px) della finestra di elaborazione orizzontale (HPW) per ottenere uno 
schema orizzontale di uscite di operatori di Sobel orizzontali (HS1-HS5); 

c3) considerare una finestra di elaborazione verticale (VPW) contenente il 
pixel (Px) in esame e pixel vicini (Pe,Pf,Pg,Ph) appartenenti ad una stessa 
colonna di pixel del pixel (Px) e precedenti e seguenti il pixel (Px); 

c4) applicare detto operatore di Sobel verticale (Vsob) a ciascun pixel (Pe- 
Ph.Px) della finestra di elaborazione verticale (VPW) per ottenere uno schema 
verticale di uscite di operatori di Sobel verticali (VS1-VS5). 
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10. Metodo secondo la rivendicazione 9, caratterizzato dal fatto che dette 
finestre di elaborazione orizzontale e verticale (HPW,VPW) contengono 
ciascuna cinque pixel e sono centrate in detto pixel (Px) in esame. 

11. Metodo secondo la rivendicazione 10, caratterizzato dal fatto che nella 
fase d2) e prevista la determinazione di gradi di appartenenza di ciascuna uscita 
dell'operatore di Sobel orizzontale (HS1-HS5) dello schema orizzontale ad un 
primo insieme fuzzy "Piccolo" e ad un primo insieme fuzzy "Grande", la 
valutazione di gradi di attivazione (cd-oc9) di un primo insieme di regole fuzzy 
ciascuna delle quali associata ad almeno uno schema orizzontale predefinito di 
uscite dell'operatore di Sobel orizzontale, la determinazione di gradi di 
appartenenza di ciascuna uscita dell'operatore di Sobel verticale (VS1-VS5) 
dello schema verticale ad un secondo insieme fuzzy "Piccolo" e ad un secondo 
insieme fuzzy "Grande" e la valutazione di gradi di attivazione (pi -(39) di un 
secondo insieme di regole fuzzy ciascuna associata ad almeno uno schema 
verticale predefinito di uscite dell'operatore di Sobel verticale. 

12. Metodo secondo la rivendicazione 11, caratterizzato dal fatto che detta 
determinazione dei gradi di appartenenza delle uscite dell'operatore di Sobel 
orizzontale (HS1-HS5) dello schema orizzontale a detti primi insiemi fuzzy 
"Piccolo" e "Grande" prevede di determinare una prima ed una seconda 
funzione di appartenenza dipendenti da dette caratteristiche globali (GF), e 
detta determinazione dei gradi di appartenenza delle uscite dell'operatore di 
Sobel verticale (VS1-VS5) dello schema verticale a detti secondi insiemi fuzzy 
"Piccolo" e "Grande" prevede di determinare una terza ed una quarta funzione 
di appartenenza dipendenti da dette caratteristiche globali (GF). 

13. Metodo secondo la rivendicazione 12, caratterizzato dal fatto che detti 
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gruppi di fattori peso predefiniti (cij) comprendono gruppi di fattori peso 
predefiniti orizzontali e gruppi di fattori peso predefiniti verticali, ciascuna 
regola fuzzy di detto primo insieme essendo associata ad un rispettivo gruppo 
di fattori peso orizzontali predefiniti, e ciascuna regola fuzzy di detto secondo 
insieme essendo associata ad un rispettivo gruppo di detti fattori peso 
predefiniti verticali. 

14. Metodo secondo la rivendicazione 13, caratterizzato dal fatto che detti 
fattori peso (H1-H5,V1-V5) comprendono fattori peso orizzontali (H1-H5) e 
fattori peso verticali (VI -V5), detti fattori peso orizzontali (H1-H5) essendo 
determinati mediante una media pesata dei fattori peso predefiniti orizzontali 
con coefficienti peso costituiti dai gradi di attivazione (al-oc9) delle regole 
fuzzy del primo insieme, e detti fattori peso verticali (VI -V5) essendo 
determinati mediante una media pesata dei fattori peso predefiniti verticali con 
coefficienti di peso costituiti dai gradi di attivazione (P1-P9) delle regole fuzzy 
del secondo insieme. 

15. Metodo secondo la rivendicazione 14, caratterizzato dal fatto che il 
valore del pixel (Px) in esame viene modificato applicando i fattori peso 
orizzontali (H1-H5) ai valori dei pixel (Pa-Pd,Px) della fmestra di elaborazione 
orizzontale (HPW) ed applicando i fattori peso verticali (VI -V5) ai valori dei 
pixel (Pe-Ph,Px) della finestra di elaborazione verticale (VPW). 

16. Dispositivo di post-elaborazione per ridurre artefatti in immagini 
digitali codificate a blocchi, caratterizzato dal fatto di comprendere: 

- primi mezzi (1) alimentati da un'immagine di ingresso (In) per stimare 
caratteristiche globali (MOUT) di un blocco di immagine (IB) in esame, dette 
caratteristiche globali fomendo informazioni su un contenuto medio di bordi di 
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immagine lungo la direzione orizzontale e verticale di detto blocco di 
immagine (IB); 

- secondi mezzi (7,8,1 1,12,D) alimentati da detta immagine di ingresso 
(In) per stimare caratteristiche locali (VS1-VS5,HS1-HS5) per ciascun pixel 
(Px) del blocco di immagine (IB) in esame, detta caratteristiche locali fomendo 
informazioni sul contenuto di bordi di immagine lungo le direzioni orizzontale 
e verticale di un'area di immagine nell'intomo di detto pixel (Px); 

- terzi mezzi (9,D,10) alimentati da dette caratteristiche globali (MOUT) e 
da dette caratteristiche locali (VS1-VS5,HS1-HS5) per modificare il valore di 
detto pixel (Px) in funzione sia di dette caratteristiche globali (MOUT) che di 
dette caratteristiche locali (VS1-VS5,HS1-HS5). 

17. Dispositivo secondo la rivendicazione 16, caratterizzato dal fatto che 
detti primi mezzi (1) comprendono mezzi (3,4) per valutare uscite di operatori 
di Sobel verticale ed orizzontale (Hsob,Vsob) applicati ai pixel di un sotto- 
blocco di immagine (ISB) intemo a detto blocco di immagine (IB) in esame, e 
mezzi accumulator! (5) per accumulare le uscite dell'operatore di Sobel 
orizzontale (Hsob) e le uscite dell'operatore di Sobel verticale (Vsob) per 
ciascun pixel del sotto-blocco di immagine (ISB). 

18. Dispositivo secondo la rivendicazione 17, caratterizzato dal fatto che 
detti secondi mezzi (7,8,1 1,12,D) comprendono quarti mezzi (7,8) per valutare 
uscite dell'operatore di Sobel verticale (VS1-VS5) applicato a detto pixel (Px) 
e a pixel verticalmente vicini (Pe,Pf,Pg,Ph) che precedono e seguono il pixel 
(Px) nella direzione verticale, e quinti mezzi (11,12,D) per valutare' le uscite 
dell'operatore di Sobel orizzontale (HS1-HS5) applicato a detto pixel (Px) e a 
pixel orizzontalmente vicini (Pa,Pb,Pc,Pd) che precedono e seguono il pixel 
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(Px) nella direzione orizzontale. 

19. Dispositivo secondo la rivendicazione 18, caratterizzato dal fatto che 
detti terzi mezzi (9,D,10) comprendono mezzi di filtraggio verticale (9,D) 
alimentati da dette uscite dell'operatore di Sobel verticale (VS1-VS5) e da 
dette caratteristiche globali (MOUT) per calcolare un valore filtrato 
verticalmente (Ex.) di detto pixel (Px) che dipende da dette caratteristiche 
globali (MOUT), da dette uscite dell'operatore di Sobel verticale (VS1-VS5) e 
dai valori di detti pixel verticalmente vicini (Pe,Pf,Pg,Ph), e mezzi di filtraggio 
orizzontale (10) alimentati da dette caratteristiche globali (MOUT), da dette 
uscite dell'operatore di Sobel orizzontale (HS1-HS5) e da detto valore filtrato 
verticalmente (Px) di detto pixel (Px) per calcolare un valore filtrato 
orizzontalmente di detto pixel (Px) che dipende da dette caratteristiche globali 
(GF), da dette uscite dell'operatore di Sobel orizzontale (HS1-HS5) e dai valori 
filtrati verticalmente (Pa, Pb , Px ,gg,Pd) di detto pixel (Px) e dei pixel 
orizzontalmente vicini (Pa,Pb,Pc,Pd). 

20. Dispositivo secondo la rivendicazione 19, caratterizzato dal fatto che 
detti mezzi di filtraggio verticale (9,D) comprendono primi mezzi di 
elaborazione fuzzy (13) alimentati da dette uscite dell'operatore di Sobel 
verticale (VS1-VS5) e da dette caratteristiche globali (MOUT) per valutare 
gradi di coincidenza (pi-(39) dello schema di uscite dell'operatore di Sobel 
verticale (VS1-VS5) con un insieme di schemi predefiniti di uscite 
dell'operatore di Sobel verticale, detti gradi di coincidenza ((3l-p9) dipendendo 
da dette caratteristiche globali (MOUT), una tabella di riferimento di 
coefficienti di filtro predefiniti (F1-F9) ciascuno associato ad almeno un 
rispettivo schema di uscite dell'operatore di Sobel verticale, mezzi di 
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composizione di coefficienti di filtro verticale (14) alimentati da detti 
coefficient! di filtro verticali predefinti per generare un insieme di coefficienti 
di filtro verticale (VI -V5) che sono una media pesata dei coefficienti di filtro 
verticale predefiniti con coefficienti di peso formati da detti gradi di 
coincidenza ((31 -f39), ed un filtro verticale (15) alimentato da detti coefficienti 
di filtro verticale (V1-V5) e dai valori del pixel (Px) e dei pixel verticalmente 
vicini (Pe,Pf,Pg,Ph) per fornire ad una uscita (FOUT) detto valore filtrato 
verticalmente del pixel (Px). 

21. Dispositivo secondo la rivendicazione 20, caratterizzato dal fatto che 
detti mezzi di filtraggio orizzontali (10) comprendono secondi mezzi di 
elaborazione fuzzy (13) alimentati da dette uscite dell'operatore di Sobel 
orizzontale (HS1-HS5) e da dette caratteristiche globali (MOUT) per valutare 
gradi di coincidenza (al-cc9) dello schema di uscite dell'operatore di Sobel 
orizzontale (HS1-HS5) con un insieme di schemi predefiniti di uscite 
dell'operatore di Sobel orizzontale, detti gradi di coincidenza (ccl-cc9) 
dipendendo da dette caratteristiche globali (MOUT), una tabella di riferimento 
di coefficienti di filtro orizzontale predefiniti (F1-F9) ciascuno associato ad 
almeno un rispettivo schema predefinito di uscite dell'operatore di Sobel 
orizzontale, mezzi di composizione di coefficienti di filtro orizzontale (14) 
alimentati da detti coefficienti di filtro orizzontale predefinti (F1-F9) per 
generare un insieme di coefficienti di filtro orizzontale (H1-H5) che sono una 
media pesata dei coefficienti di filtro orizzontale predefiniti con coefficienti 
peso costituiti da detti gradi di coincidenza (al-cc9), ed un filtro orizzontale 
(15) alimentato da detti coefficienti di filtro orizzontale (H1-H5) e dei valori 
filtrati verticalmente (Pa.Pb,Px,Pc,Pd) del pixel (Px) e dei pixel 
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orizzontalmente vicini (Pa,Pb,Pc,Pd) per fomire ad una uscita (FOUT) detto 
valore filtrato orizzontalmente (Out) del pixel (Px). 
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"Metodo di post-elaborazione per ridurre artefatti in irrtmagini digitali 
codificate a blocchi, e dispositivi di post-elaborazione per attuare tale metodo." 

$ )j< $ j(e s(e 

RIASSUNTO 

Un metodo di post-elaborazione per ridurre gli artefatti in immagini 
digitali codificate a blocchi, caratterizzato dal fatto di prevedere: a) la 
suddivisione di un immagine in ingresso (In) in una pluralita di blocchi di 
immagine (IB); b) per ciascun blocco di immagine (IB), la stima di 
caratteristiche globali (GF) di detto blocco di immagine che forniscono 
informazioni sul contenuto medio di bordi di immagine lungo le direzioni 
orizzontale e verticale di detto blocco di immagine (IB); c) per ciascun pixel 
(Px) di un blocco di immagine (IB) in esame, la stima di caratteristiche locali 
(LF) relative a detto pixel (Px) che forniscono informazioni sul contenuto di 
bordi di immagine lungo la direzione orizzontale e verticale di un'area di 
immagine nell'mtorno di detto pixel (Px); d) la modifica del valore di detto 
pixel (Px) in funzione sia di dette caratteristiche globali (GF) del blocco di 
immagine (IB) cui appartiene detto pixel (Px), sia delle caratteristiche locali 
(LF) dell'area di immagine nell'intorno di detto pixel (Px). (Figura 1). 



1 



1 



"Post-processing method for reducing artifacts in block-coded digital 
images, and post-processing device for actuating such method." 

* jfc 3fc * 

DESCRIPTION 

5 The present invention relates digital image coding. More precisely, the 

invention relates to a post-processing method for reducing artifacts in block- 
coded digital images, and to a post-processing device suitable for actuating 
such a method. 

With the diffusion of digital communication systems, digital images are 
10 more and more used. This has led to the diffusion of still and video cameras 

with digital acquisition and processing capability. 

In order to better exploit storage means and transmission bandwidth, 
digital image compression standards have been developed, such as JPEG for 
still images, and MPEG-1 and MPEG-2 for digital television image 
15 sequences. 

The above-referred compression standards provides for block-coding 
based on Discrete Cosine Transform (DCT). A digital image is divided into 
blocks of pixels, and each block is encoded independently from the others. 
DCT coefficients for the pixels of each block are evaluated and a 
20 quantization matrix is applied to the DCT coefficients to reduce the 

information to be stored or transmitted. When the image is to be displayed, it 
must be decoded in advance. 

Due to the quantization process, these image compression methods are 
lossy, i.e. they cause a loss of information in the decoded image with respect 
25 to the original image. The decoded image can thus present noticeable 

degradation, mainly consisting of two kinds of artifacts known in the art 
under the names of "grid noise" and "staircase noise". 

In order to reduce the image degradation, post-processing methods of 
the decoded image have been proposed which allow for attenuating grid 
30 noise and staircase noise. 
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In view of the state of the art described, it is an object of the present 
invention to provide a new post-processing method for reducing artifacts in 
block-coded digital images. 

According to the present invention, such object is attained by means of a 
5 post-processing method for reducing artifacts in block-coded digital images, 

characterized by providing for: 

a) dividing an input image into a plurality of image blocks; 

b) for each image block, estimating global features of said image 
block providing information on an average content of image 

1 o edges along the horizontal and vertical directions of said image 

block; 

c) for each pixel of an image block under examination, estimating 
local features for said pixel providing information on the content 
of image edges along the horizontal and vertical directions of an 

15 image area around said pixel; 

d) modifying the value of said pixel according to both said global 
features of the image block to which said pixel belongs and said 
local features of the image area around said pixel. 

Also according to the invention, there is also provided a post-processing 
20 device for actuating said method, characterized by comprising: 

first means supplied with an input image for estimating global 
features of an image block under examination, said global features 
providing information on an average content of image edges along 
the horizontal and vertical directions of said image block; 
25 - second means supplied with said input image for estimating local 

features for each pixel of the image block under examination, said 
local features providing information on the content of image edges 
along the horizontal and vertical directions of an image area around 
said pixel; 

30 - third means supplied with said global features and said local features 
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for modifying the value of said pixel according to both said global 
features and said local features. 
The features and advantages of the present invention will be made 
apparent from the following detailed description of an embodiment thereof, 
5 illustrated as a non-limiting example in the annexed drawings, wherein: 

Figure 1 is a schematic block diagram illustrating the principle of 
operation of a method according to the present invention; 

Figure 2 shows a digital image divided into image blocks; 
Figure 3 shows in detail an image block into which the digital image of 
1 0 Figure 2 is divided; 

Figure 4 shows an array of pixels of the image block of Figure 3; 
Figure 5 shows an image sub-block of the image block of Figure 3 used 
for evaluating global features of the image block; 

Figure 6 shows an horizontal processing window used for evaluating 
15 local features in the horizontal direction for a generic pixel of the image 

block; 

Figure 7 shows a vertical processing window used for evaluating local 
features in the vertical direction for said generic pixel; 

Figures 8 and 9 shows two membership functions used to perform a 
20 fuzzy computation; 

Figure 10 is a block diagram of a device according to the present 
invention; 

Figure 1 1 shows the structure of two blocks of the device of Figure 10; 

and 

25 Figure 12 is a block diagram of other two blocks of the device of Figure 

10. 

With reference to Figure 1, there is shown a block diagram illustrating 
the principle of operation of the post-processing method according to the 
present invention. An input decoded compressed digital image In is supplied 
30 to a Feature Extraction block FE. Block FE provides for analyzing the image 
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to evaluate global and local features thereof. The global and local features, 
respectively GF and LF, of the image In evaluated by block FE are supplied 
as inputs to a Fuzzy Process block FUZZY which, according to fuzzy rules, 
determines parameters FA suitable for determining the kind of filtering to be 
5 performed, in accordance to the global and local features GF and FL of the 

image. The parameters FA calculated by block FUZZY are supplied to a 
Filter Composition block FC which according to said parameters FA 
determines the type of filtering to performed out of a set of predefined filters 
(block FS). Filter parameters FP determined" by block FC are then supplied 

10 to a Processing block PROC, also supplied directly with the input image In, 

which performs the filtering of the input image In according to the filter 
parameters FP to provide a post-processed output image Out. 

It appears that the kind of filtering to be performed on the decoded input 
image In is chosen after an estimation of the global and local features of the 

15 decoded input image. For image areas near grid noise and near an edge, a 

low-pass filtering is performed, to reduce both staircase noise and grid noise. 
For areas containing fine details (image edges and texture), no filtering is 
performed. Thus, the method according to the present invention provides for 
performing a non-linear adaptive filtering on the pixels of the decoded 

20 image. 

The principle of operation outlined above will be now described in 
detail. 

As shown in Figure 2, the input image In is partitioned into image 
blocks IB, each containing an equal number of pixels. A typical dimension 
25 of the blocks is 8*8 pixels (Figure 3), but this is not however to be intended 

as a limitation, since other block dimensions are suitable. 

The image blocks IB of the input image In are scanned line by line 
starting from the top-left block to the bottom-right one. For each image 
block IB, the Feature Extraction block FE in Figure 1 determines the global 
30 and local features GF and LF. 
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Global features of the image block IB under examination are determined 
by applying horizontal and vertical Sobel operators: 



horizontal Sobel operator (Hsob): 



vertical Sobel operator (Vsob) 



hll hl2 hl3 
h21 h22 h33 
h31 h32 h33 

vll vl2 vl3^ 
v21 v22 v23 
v31 v32 v33 
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to each pixel belonging to an image sub-block internal to the image block 
IB. For example, the following Sobel operators: 
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are applied to a 6*6 pixel image sub-block ISB (gray area in Figure 5). As 
shown in Figure 4, for each pixel Pi of the image sub-block ISB a 3*3 array 
of neighboring pixels M centered in pixel Pi is considered, and the values of 
the pixels of said array M are multiplied by the coefficients of the horizontal 
and vertical Sobel operators, to obtain: 

Hsob = (P3+P8+2*P5) - (P1+P6+2*P4), 

Vsob = (P6+P8+2*P7) - (P1+P3+2*P2), 
where P1-P4 and P5-P8 are the values of the pixels (gray levels). 

As known, horizontal and vertical Sobel operators perform a filtering 
capable of detecting edges along the horizontal and vertical direction, 
respectively. 

The output values of the horizontal Sobel operators calculated for the 
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pixels of image sub-block ISB are accumulated to obtain an accumulated 
value Acc(Hsob), and the output values of the vertical Sobel operators 
calculated for the pixels of image sub-block ISB are accumulated to obtain 
an accumulated value Acc(Vsob). Acc(Hsob) gives the high-frequency 
5 content in the horizontal direction (vertical edges) of the image block IB. 

Acc(Vsob) gives the high-frequency content in the vertical direction 
(horizontal edges) of image block IB. Thus, Acc(Hsob) and Acc(Vsob) 
respectively provide the degree of "edgeness" of the image block under 
examination in the vertical and horizontal directions. It is to be noted that in 
1 0 order to evaluate the global features GF of the image block IB under 

examination, only the pixels belonging to this block are considered (by 
applying 3*3 Sobel operators to the 6*6 image sub-block ISB, it is not 
necessary to take into consideration pixels belonging to neighboring image 
blocks). 

15 Global features GF of the image block under examination can be formed 

by the accumulated values Acc(Hsob) and Acc(Vsob). Alternatively, the 
global features GF of the image block can be formed by an average of the 
accumulated values Acc(Hsob) and Acc(Vsob), to obtain the average 
number of edges in the horizontal and vertical directions in the image block 

20 under examination. 

Local features LF of the image block IB are estimated according to the 
following method. All the pixels of the image block IB under examination 
are scanned line by line starting from the top-left pixel down to the bottom- 
right one. To each pixel a horizontal processing window is applied: a 

25 prescribed number of pixels respectively preceding and following the pixel 

under consideration and belonging to the same image line of the pixel under 
consideration is considered. A suitable horizontal processing window FfPW 
is shown in Figure 6, which is a horizontal 1*5 processing window: for a 
given pixel, the two preceding pixels Pa, Pb and the two following pixels Pc, 

30 Pd belonging to the same line are considered. In Figure 6 there is shown by 
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way of example the horizontal processing window HPW associated to the 
first pixel Px of the image block. It should be noted that not only the pixels 
of the image block IB under examination are considered, but also pixels 
belonging to neighboring image blocks; this is for example the case of the 
first, second, seventh and eight pixel of each line of pixels of the image 
block IB under examination. 

The horizontal Sobel operator Hsob previously mentioned is applied to 
each pixel Pa, Pb, Px, Pc, Pd in the horizontal processing window HPW, to 
obtain five output values HS1-HS5. Values HS1-HS5 provide the local 
features in the horizontal direction for the pixel under examination Px, i.e. 
the high-frequency content in the horizontal direction of the image region 
around the pixel under examination. 

Similarly, a vertical processing window is applied to each pixel of the 
image block IB. The vertical processing window is formed by the pixel 
under consideration Px, and a prescribed number of pixels belonging to the 
same column as and preceding and following the pixel under consideration; 
for example, as shown in Figure 7 the vertical processing window VPW can 
have dimensions identical to the horizontal processing window HPW (5*1), 
and thus contains two pixels Pe, Pf preceding pixel Px and two pixels Pg, Ph 
following pixel Px in the vertical direction. 

The vertical Sobel operator Vsob previously mentioned is then applied 
to each pixel of the vertical processing window VPW to obtain five output 
values VS1-VS5. Values VS1-VS5 form the local features in the vertical 
direction for the pixel under examination, i.e. the high-frequency content in 
the vertical direction of an image region around the pixel under examination. 

The global features GF for the image block IB under examination (i.e., 
the two accumulated values Acc(Hsob) and Acc(Vsob) or, in alternative, the 
average value of Acc(Hsob) and Acc(Vsob)) and the local features LF for 
the pixel under examination inside said image block (the ten values HS 1- 
HS5 and VS1-VS5) are then supplied to the Fuzzy Process block FUZZY. 
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This block provides for evaluating the degree of membership of a generic 
value HSi and Vsi (i=1..5) to two fuzzy sets "Small" and "Big"; these degree 
of membership can be evaluated by applying to HSi, VSi the membership 
functions depicted in Figures 8 and 9. In these figures, Thl and Th2 are 
values depending on the global features GF of the image block under 
examination, i.e. on the accumulated values Acc(Hsob) and Acc(Vsob) or on 
the average of the accumulated values; in the first case, Thl and Th2 are 
different for the Hsi and Vsi values; in the second case, Thl and Th2 are the 
same for Hsi and Vsi values. 

Fuzzy rules having as antecedents the degrees of membership of the 
output values HSi and VSi to the two fuzzy sets "Small" and "Big" are then 
evaluated. This means that 32 rules are to be evaluated for both the 
horizontal and vertical directions. However, all those fuzzy rules having the 
same consequence are synthesized in one rule only by an else operator. In 
this way, the system complexity is reduced, and a total of nine rules for each 
direction have to be evaluated. 

The following fuzzy rules are applied to the five values HS1-HS5 
associated to the horizontal direction: 

1. If HSI is Small and HS2 is Small and HS3 is Small and HS4 is Small 
and HS5 is Small, then al is Big; 

2. If HSI is Small and HS2 is Small and HS3 is Small and HS4 is Small 
and HS5 is Big, then ct2 is Big; 

3. If HSI is Small and HS2 is Small and HS3 is Small and HS4 is Big 
and HS5 is Small, then cc3 is Big; 

4. If HSI is Small and HS2 is Small and HS3 is Small and HS4 is Big 
and HS5 is Big, then cc4 is Big; 

5. If HSI is Small and HS2 is Big and HS3 is Small and HS4 is Small 
and HS5 is Small, then ct5 is Big; 

6. If HSI is Big andHS2 is Small and HS3 is Small and HS4 is Small 
and HS5 is Small, then a.6 is Big; 
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7. IfHS 1 is Big and HS2 is Small and HS3 is Small arc J HS4 is Small 
a«J //S5 is Big, f/ie>7 a7 z'5 Big; 

8. IfHSl is Big andHS2 is Big and HS3 is Small and HS4 is Small and 
HS5 w Small, zTze/i a8 is Big. 

5 The activation level of each rule depends on the degree of memberships 

of the pattern of output values HSi of the horizontal Sobel operator applied 
to the five pixels of the horizontal processing window HPW; the degree of 
memberships depend in turn on the global features GF of the image block to 
which the pixel under examination belongs. The activation level of the else 
10 (ninth) rule is computed as ccelse = (1 - ocave), where aave is the average 

activation degree of fuzzy rules 1 to 8. al to a.8 and aelse, and a similar set 
of nine activation degrees for the fuzzy rules applied to values VS1-VS5) 
form the output FA of the fuzzy process block FUZZY in Figure 1 . 

Each one of the above-listed rules is associated to a respective set of 
15 predefined filter parameters, which are stored as a look-up table in block FS 

of Figure 1. Suitable predefined filter parameter sets are for example: 
Rule 1: (cl 1=1.0, cl2=1.0, cl3=1.0, cl4=1.0, cl5=1.0) if the pixel under 
examination lies outside the image sub-block ISB, and 
(cl 1=0.0', cl2=1.0, cl3=1.0, cl4=1.0, cl5=0.0) if the pixel under 
20 examination lies inside the image sub-block ISB; 

Rule 2: (c21=0.5, c22=1.0, c23=1.0, c24=1.0, c25=0.0); 
Rule 3: (c31=0.5, c32=1.0, c33=1.0, c34=0.0, c35=0.0); 
Rule 4: (c41=0.5, c42=1.0, c43=1.0, c44=0.0, c45=0.0); 
Rule 5: (c5 1=0.0, c52=0.0, c53=1.0, c54=1.0, c55=0.5); 
25 Rule 6: (c61=0.0, c62=1.0, c63=1.0, c64=1.0, c65=0.5); 

Rule 7: (c71=0.0, c72=1.0, c73=1.0, c74=1.0, c75=0.0); 
Rule 8: (c81=0.0, c82=0.0, c83=1.0, c84=1.0, c85=0.5); 
. Else rule: (c91=0.0, c92=0.0, c93=1.0, c94=0.0, c95=0.0). 

The parameters FP of the filter to be applied to the five pixels of the 
30 horizontal processing window HPW are calculated as a weighted average of 
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the nine filters described above, with weight factors formed by the activation 
degrees al to a8 and aelse of the respective fuzzy rules. 

Assuming that ai is the activation degree of the i-th fuzzy rule (i=1..9), 
the ninth fuzzy rule being the else fuzzy rule (a9 = aelse), and cij are the 
coefficients of the i-th filter (i=1..9, j=1..5), the weight factor applied to the 
i-th filter, associated to the i-th fuzzy rule is: 

Fi = ai • cij 

and the coefficients Hj of the final horizontal filter to be applied to the pixels 

of the horizontal processing window HPW are given by: 
9 

Hj = & : 

N 

where N is a normalization factor. 

The horizontally- filtered value Px of the pixel Px under examination (at 
the center of the horizontal processing window) is then calculated as a 
weighted average of the values of the pixels Pa, Pb, Px, Pc and Pd belonging 
to the horizontal processing window HPW, with weight factors formed by 
the coefficients Hj : 

Px = Hl*Pa+ H2*Pb+ H3*Px+ H4*Pc+ H5*Pd. 

Similar calculations are performed for the vertical direction, starting 
from the output values VS1-VS5 of the vertical Sobel operators applied to 
the pixels Pe, Pf, Px, Pg and Ph in the vertical processing window VPW. 
The coefficients Vj (j=1..5) of the filter for the vertical direction are 
calculated in a way completely similar to that used for determining the 
coefficients Hj: 
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where pi (i=1..9) are the activation degrees of nine fuzzy rules for the 
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vertical direction (similar to those listed above for the horizontal direction) 
and cij (j=1..5) now are the predefined filter parameters associated to the i-th 
fuzzy rule for the vertical direction. The coefficients Vj are then applied to 
the pixels in the vertical processing window VPW to calculate a weighted 
average of the same. The filtered value of the pixel Px under examination, 
filtered in both the horizontal and vertical direction, is provided at the output 
Out of the processing block PROC. 

The value of the pixel Px under examination to be multiplied by the 
vertical filter coefficient V3 can be the value Px obtained after having 
applied to the pixels in the horizontal processing window HPW the 
horizontal filter Hj (j=1..5): 

Out = Vl*Pe + V2*Pf + V3*Px + V4*Pg + V5*Ph. 

Alternatively, it is possible to evaluate first the vertically- filtered value 
Px of the pixel under examination: 

Ex = Vl*Pe + V2*Pf + V3*Px + V4*Pg + V5*Ph, 
and then performing the filtering in the horizontal direction applying to this 
value the respective coefficient H3 of the horizontal filter Hj : 
Out = HI *Pa + H2*Pb + H3*Ex + H4*Pc + H5*Pd. 

The sequence is of no importance, the important thing to be underlined being 
that at the end of the process the value of the pixel under examination is the 
result of both an horizontal and a vertical filtering. 

Figure 10 is a block diagram of a device suitable for actuating the 
method previously described. The device comprises two main blocks: a first 
block 1 evaluates the global features GF of the image blocks IB the image to 
be post-processed is divided in, and a second block 2 evaluates the local 
features LF of the pixels of the image and performs the filtering according to 
both the global features and the local features. 

It is assumed that the image to be post-processed is scanned line by line 
in a sequential order. Signal In is a stream of pixels of the input image 
scanned line by line. Block 1 is supplied with signal In; signal In also 
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supplies a cascade of two line memories LM1 and LM2 whose outputs 
supplies block 1. 

Inside block 1, signal In and the outputs of line memories LM1 and 
LM2 supply a block 3 of pixel delays suitable for implementing a 3*3 pixel 
window which is used to calculate horizontal and vertical Sobel operators 
for the pixels of the 6*6 image sub-block ISB inside each image block IB. 
Block 3 supplies a block 4 which calculates the outputs Hsob and Vsob of 
the horizontal and vertical Sobel operators for those pixels of the current 
image line belonging to the 6*6 image sub-blocks ISB of each image block 
IB. The outputs Hsob and Vsob of block 4 are supplied to an accumulator 
block 5 wherein they are accumulated. After eight image lines, i.e. a line of 
image blocks IB, have been scanned, the accumulated values Acc(Hsob), 
Acc(Vsob) (or alternatively the average thereof) for each image block IB are 
stored in a memory block 6. 

The output of line memory LM2 supplies a cascade of eight further line 
memories LM3-LM10. Block 2 is supplied in parallel with the outputs of 
line memories LM4-LM10. In this way, evaluation of the local features and 
calculation of the filter parameters starts after block 1 has estimated the 
global features GF for a line of image blocks IB. 

Inside block 2, a block 7 of pixel delays is supplied with the outputs of 
line memories LM4-LM10; by means of the line memories LM4-LM10 and 
the pixel delay block 7 it is possibile to implement the 5*1 vertical 
processing window VPW. The outputs L4-L10 of the pixel delays block 7 
supply a block 8 which applies the vertical Sobel operator to each pixel 
inside the vertical processing window VPW. To avoid the use of further line 
memories, a parallel approach is preferred providing for calculating five 
vertical Sobel operators in parallel; the outputs of the five vertical Sobel 
operators VS1-VS5 are supplied to a vertical fuzzy filter 9, which is also 
supplied with the outputs L6-L10 of the pixel delay block 7 and the output 
MOUT of the memory MEM of block 1. MOUT supplies the global features 
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GF of the image block IB currently processed by block 2, i.e. the 
accumulated value Acc(Vsob) or, alternatively, the average of Acc(Vsob) 
and Acc(Hsob). The vertical fuzzy filter block 9 evaluates the degree of 
membership of values VS1-VS5 to the fuzzy sets "Small" and "Big" taking 
5 into account the global features provided by MOUT, evaluates the activation 

levels of the nine fuzzy rules for the vertical direction, calculates the 
coefficients Vj (j=1..5) of the vertical filter and applies the vertical filter 
coefficients Vj to the five pixels Pe, Pf, Px, Pg, Ph in the vertical processing 
window VPW, to calculate the vertically-filtered value Ex. of the pixel in the 

10 middle of the vertical processing window. The output of the vertical fuzzy 

filter block 9 forms the vertically- filtered value Ed. of pixel Pd in the 
horizontal processing window HPW shown in Figure 6 and supplies directly 
a horizontal fuzzy filter block 10; the output Ed of block 9 also supplies a 
cascade of four pixel delays D whose outputs respectively form the 

1 5 vertically- filtered values Eg, 2x, Eh, Ej| of the pixels Pc, Px, Pb, Pa in the 

horizontal processing window HPW and supply the horizontal fuzzy filter 
block 10. 

In parallel to blocks 8 and 9, the outputs L7-L9 of the pixel delay block 
7 supply a block 1 1 which applies the horizontal Sobel operators to the 

20 pixels inside the horizontal processing window HPW. Differently from the 

vertical sobel operators, only one horizontal sobel operator is calculated at a 
time; a compensation delay block 12 introduces a delay for compensating 
the processing delay of the vertical fuzzy filter block 9. The output of block 
12, forming the output of the horizontal Sobel operator HS5 applied to pixel 

25 Pd of the horizontal processing window in Figure 6, supplies the horizontal 

fuzzy filter block 10 and a cascade of four pixel delays D, the outputs 
thereof forming the values HS4, HS3, HS2 and HS1 and suppling the 
horizontal fuzzy filter block 10. The horizontal fuzzy filter block 10, which 
is also supplied by the output MOUT of the memory block MEM in block 1 

30 providing the value Acc(Hsob) (or alternatively the average of values* 
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Acc(Hsob) and Acc(Vsob)), evaluates the degree of membership of values 
HS1-HS5 to the fuzzy sets "Small" and "Big" according to the value of the 
global features GF provided by MOUT, evaluates the activation levels of the 
nine fuzzy rules described above for the filtering in the horizontal direction, 
calculates the coefficients Hj of the horizontal filter and applies the 
parameters Hj to the vertically-filtered values Ej, Eb_, Ex., Ee, Ed of the pixels 
Pa, Pb, Px, Pc, Pd in the horizontal processing window HPW to obtain the 
horizontally- and vertically-filtered value Out of the pixel Px under 
examination. 

A control circuit CTRL controls the operation of blocks 1, 2 and the line 
memories LM1-LM10. 

Figure 1 1 shows the structure of the vertical and horizontal Sobel 
operator blocks 8 and 1 1 of Figure 10. They are composed in a 
straightforward way by adders. 

Figure 12 shows the structure of both the vertical fuzzy filter block 9 
and the horizontal fuzzy filter block 10. XI -X5 are the vertical or, 
respectively, horizontal Sobel operator outputs VS1-VS5 and HS1-HS5. XI- 
X5 are supplied to a fuzzy rule evaluation block 1 3 which evaluates the 
activation degrees [31 -[39 of the nine fuzzy rules for the vertical direction or, 
respectively, the activation degrees ccl-cc9 of the nine fuzzy rules for the 
horizontal direction. The activation degrees evaluated by block 1 3 are 
supplied to a look-up table of respective predefined filter parameters F1-F9 
(forming block FS in Figure 1), and the outputs of the look-up table, i.e. the 
predefined filter parameters cij multiplied by the activation degree of the 
respective fuzzy rule, are supplied to a filter composition block 14 which 
calculates the coefficients V1-V5 or, respectively, H1-H5, of the vertical or, 
respectively, horizontal filter. Said coefficients are then supplied to a 
processing block 15 which is also supplied with the pixel values PXS (L6- 
L10 or, respectively, Es, Eh, Ex, Ee, Ed in Figure 10); block 15 applies the 
filter coefficients to the pixel values to obtain the filtered value of the pixel 



under examination Px. 
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CLAIMS 

1 . Post-processing method for reducing artifacts in block-coded digital 
images, characterized by providing for: 

a) dividing an input image (In) into a plurality of image blocks (IB); 
5 b) for each image block (IB), estimating global features (GF) of said 

image block providing information on an average content of image edges 
along the horizontal and vertical directions of said image block (IB); 

c) for each pixel (Px) of an image block (IB) under examination, 
estimating local features (LF) for said pixel (Px) providing information on 

1 0 the content of image edges along the horizontal and vertical directions of an 

image area near said pixel (Px); 

d) modifying the value of said pixel (Px) according to both said 
global features (GF) of the image block (IB) to which said pixel (Px) 
belongs and said local features (LF) of the image area near said pixel (Px). 

15 2. Method according to claim 1, characterized in that step d) provides 

for: 

dl) defining a set of predefined local features; 

d2) determining degrees of coincidence (al-a9,(31-(39)of said local 
features (LF) of the image area around said pixel (Px) with each 
20 predefined local features of said set, said degrees of coincidence (ccl- 

a9,pi-P°) depending on said global features (GF) of the image block 
(DB) to which said pixel (Px) belongs; 

d3) making the value of said pixel (Px) equal to a weighted average 
( Px , Px ) of the value of said pixel (Px) and of the values of 
25 neighboring pixels (Pa-Pd,Pe-Pf), with weight factors (H1-H5.V1- 

V5) depending on said degrees of coincidence of said local features 
(LF) with each of said predefined local features. 
3. Method according to claim 2, characterized in that said determining 
the degrees of coincidence (al-a9,(31-P9) in step d2) provides for 
30 performing a fuzzy calculation. 
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4. Method according to claim 3, characterized in that each of said 
predefined local features is associated to a respective group of predefined 
weight factors (cij), and each of said weight factors (H1-H5,V1-V5) is 
calculated as a weighted average of corresponding predefined weight factors 
(cij) of said groups with weight coefficients (ctl-ct9,pi-P9) depending on 
said degrees of coincidence of said local features (LF) with each of said 
predefined local features. 

5. Method according to claim 4, characterized in that said estimating 
global features (GF) of the image block (IB) under examination provides for 
applying horizontal and vertical Sobel operators (Hsob,Vsob) to pixels 
belonging to an image sub-block (ISB) internal to said image block (IB) 
under examination. 

6. Method according to claim 5, characterized in that said estimating 
global features (GF) of the image block (IB) under examination provides for 
adding outputs of the horizontal Sobel operators (Hsob) applied to each pixel 
of said image sub-block (ISB) to obtain an accumulated output of horizontal 
Sobel operators (Acc(Hsob)), and adding outputs of the vertical Sobel 
operators (Vsob) applied to each pixel of said image sub-block (ISB) to 
obtain an accumulated output of vertical Sobel operators (Acc(Vsob)). 

7. Method according to claim 6, characterized in that said global 
features (GF) of the image block (IB) under examination are formed by said 
accumulated outputs of the horizontal and vertical Sobel operators 

( Acc(Hsob), Acc(Vsob)) . 

8. Method according to claim 6, characterized in that said global 
features (GF) of the image block (IB) under examination are formed by an 
average of said accumulated outputs of the horizontal and vertical Sobel 
operators (Acc(Hsob),Acc(Vsob)). 

9. Method according to claim 7 or 8, characterized in that said 
estimating local features (LF) for a pixel (Px) of the image block (IB) under 
examination provides for: 



cl) considering an horizontal processing window (HPW) containing 
the pixel (Px) under examination and neighboring pixels (Pa,Pb,Pc,Pd) 
belonging to a same image line as the pixel (Px) and preceding and 
following the pixel (Px); 

c2) applying said horizontal Sobel operator (Hsob) to each pixel (Pa- 
Pd,Px) of the horizontal processing window (HPW) to obtain a horizontal 
pattern of horizontal Sobel operator outputs (HS1-HS5); 

c3) considering a vertical processing window (VPW) containing the 
pixel (Px) under examination and neighboring pixels (Pe,Pf,Pg,Ph) 
belonging to a same column of pixels as the pixel (Px) and preceding and 
following the pixel (Px); 

c4) applying said vertical Sobel operator (Vsob) to each pixel (Pe- 
Ph,Px) of the vertical processing window (VPW) to obtain a vertical partem 
of vertical Sobel operator outputs (VS1-VS5). 

10. Method according to claim 9, characterized in that said horizontal 
and vertical processing windows (HPW, VPW) contains each one five pixels 
and are centered at said pixel (Px) under examination. 

1 1 . Method according to claim 1 0, characterized in that step d2) 
provides for determining degrees of membership of each horizontal Sobel 
operator output (HS1-HS5) of the horizontal pattern to a first fuzzy set 
"Small" and to a first fuzzy set "Big", evaluating activation degrees (ocl-a9) 
of a first set of fuzzy rules each one associated with at least one predefined 
horizontal pattern of horizontal Sobel operator outputs, determining degrees 
of membership of each vertical Sobel operator output (VS1-VS5) of the 
vertical pattern to a second fuzzy set "Small" and to a second fuzzy set 
"Big" and evaluating activation degrees ((31-P9) of a second set of fuzzy 
rules each one associated with at least one predefined vertical partem of 
vertical Sobel operator outputs. 

12. Method according to claim 11, characterized in that said 
determining degrees of membership of the horizontal Sobel operator outputs 
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(HS1-HS5) of the horizontal pattern to said first fuzzy sets "Small" and 
"Big" provides for determining a first and a second membership functions 
depending on said global features (GF), and said determining degrees of 
membership of the vertical Sobel operator outputs (VS1-VS5) of the vertical 
5 pattern to said second fuzzy sets "Small" and "Big" provides for determining 

a third and fourth membership functions depending on said global features 
(GF). 

13. Method according to claim 12, characterized in that said groups of 
predefined weight factors (cij) comprise groups of predefined horizontal 
10 weight factors and groups of predefined vertical weight factors, each fuzzy 

rule of said first set being associated with a respective one of said groups of 
predefined horizontal weight factors, and each fuzzy rule of said second set 
being associated with a respective one of said groups of predefined vertical 
weight factors. 

15 14. Method according to claim 13, characterized in that said weight 

factors (H1-H5,V1-V5) comprise horizontal weight factors (H1-H5) and 
vertical weight factors (VI -V5), said horizontal weight factors (H1-H5) 
being determined by making a weighted average of the predefined horizontal 
weight factors with weight coefficients being formed by the activation 

20 degrees (al-a9) of the fuzzy rules of the first set, and said vertical weight 

factors (VI -V5) being determined by making a weighted average of the 
predefined vertical weight factors with weight coefficients formed by the 
activation degrees ([31 -(39) of the fuzzy rules of the second set. 

15. Method according to claim 14, characterized in that the value of the 

25 pixel (Px) under examination is modified by applying the horizontal weight 

factors (H1-H5) to the values of the pixels (Pa-Pd,Px) in the horizontal 
processing window (HPW) and applying the vertical weight factors (V1-V5) 
to the values of the pixels (Pe-Ph,Px) in the vertical processing window 
(VPW). 

30 16. Post-processing device for reducing artifacts in block-coded digital 
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images, characterized by comprising: 

first means (1) supplied with an input image (In) for estimating 
global features (MOUT) of an image block (IB) under examination, 
said global features providing information on an average content of 
5 image edges along the horizontal and vertical directions of said 

image block (IB); 

second means (7,8,1 1,12,D) supplied with said input image (In) for 
estimating local features (VS1-VS5,HS1-HS5) for each pixel (Px) of 
the image block (IB) under examination, said local features 
10 providing information on the content of image edges along the 

horizontal and vertical directions of an image area around said pixel 
(Px); 

third means (9,D,10) supplied with said global features (MOUT) and 
said local features (VS1-VS5,HS1-HS5) for modifying the value of 
1 5 said pixel (Px) according to both said global features (MOUT) and 

said local features (VS1-VS5,HS1-HS5). 

17. Device according to claim 16, characterized in that said first means 
(1) comprises means (3,4) for evaluating horizontal and vertical Sobel 
operator outputs (Hsob,Vsob) of horizontal and vertical Sobel operators 

20 applied to pixels of an image sub-block (ISB) internal to said image block 

(IB) under examination, and accumulator means (5) for accumulating the 
horizontal Sobel operator outputs (Hsob) and the vertical Sobel operator 
outputs (Vsob) for each pixel of the image sub-block (ISB). 

18. Device according to claim 17, characterized in that said second 
25 means (7,8,1 1,1 2,D) comprises fourth means (7,8) for evaluating vertical 

Sobel operator outputs (VS1-VS5) of vertical Sobel operators (Vsob) 
applied to said pixel (Px) and to vertically-neighboring pixels (Pe,Pf,Pg,Ph) 
preceding and following the pixel (Px) in the vertical direction, and fifth 
means (11,12,D) for evaluating horizontal Sobel operator outputs (HS1- 
30 HS5) of horizontal Sobel operators (Hsob) applied to said pixel (Px) and to 



horizontally-neighboring pixels (Pa,Pb,Pc,Pd) preceding and following the 
pixel (Px) in the horizontal direction. 

19. Device according to claim 18, characterized in that said third means 
(9,D,10) comprises vertical filter means (9,D) supplied with said vertical 
Sobel operator outputs (VS1-VS5) and said global features (MOUT) for 
calculating a vertically-filtered value ( Px") of said pixel (Px) depending on 
said global features (MOUT), said vertical Sobel operator outputs (VS1- 
VS5) and the values of said vertically-neighboring pixels (Pe,Pf,Pg,Ph), and 
horizontal filter means (10) supplied with said global features (MOUT), said 
horizontal Sobel operator outputs (HS1-HS5) and said vertically-filtered 
value ( Px) of said pixel (Px) for calculating a horizontally-filtered value of 
said pixel (Px) depending on said global features (MOUT), said horizontal 
Sobel operator outputs (HS1-HS5) and the vertically-filtered values 
( Pa . Pb . Px .Pc T Pd ) of said pixel (Px) and said horizontally-neighboring pixels 
(Pa,Pb,Pc,Pd). 

20. Device according to claim 19, characterized in that said vertical 
filter means (9,D) comprises first fuzzy computation means (13) supplied 
with said vertical Sobel operator outputs (VS1-VS5) and said global features 
(MOUT) for evaluating degrees of coincidence ((31-p9) of the pattern of 
vertical Sobel operator outputs (VS1-VS5) with a set of predefined patterns 
of vertical Sobel operator outputs, said degrees of coincidence ((31 -(39) 
depending on said global features (MOUT), a look-up table of predefined 
vertical filter coefficients (F1-F9) each one associated to at least one 
respective predefined pattern of vertical Sobel operator outputs, a vertical 
filter coefficients composition means (14) supplied with said predefined 
vertical filter coefficients for generating a set of vertical filter coefficients 
(V1-V5) which are a weighted average of the predefined vertical filter 
coefficients with weight coefficients formed by said degrees of coincidence 
((31 -(39), and a vertical filter (15) supplied with said vertical filter 
coefficients (VI -V5) and the values of the pixel (Px) and the vertically- 



neighboring pixels (Pe,Pf,Pg,Ph) for providing at an ouput (FOUT) said 
vertically-filtered value of the pixel (Px). 

2 1 . Device according to claim 20, characterized in that said horizontal 
filter means (10) comprises second fuzzy computation means (13) supplied 
with said horizontal Sobel operator outputs (HS1-HS5) and said global 
features (MOUT) for evaluating degrees of coincidence (al-a.9) of the 
pattern of horizontal Sobel operator outputs (HS1-HS5) with a set of 
predefined patterns of horizontal Sobel operator outputs, said degrees of 
coincidence (ocl-a9) depending on said global features (MOUT), a look-up 
table of predefined horizontal filter coefficients (F1-F9) each one associated 
to at least one respective predefined pattern of horizontal Sobel operator 
outputs, a horizontal filter coefficients composition means (14) supplied 
with said predefined horizontal filter coefficients (F1-F9) for generating a set 
of horizontal filter coefficients (H1-H5) which are a weighted average of the 
predefined horizontal filter coefficients with weight coefficients formed by 
said degrees of coincidence (ccl-a9), and a horizontal filter (15) supplied 
with said horizontal filter coefficients (H1-H5) and the vertically-filtered 
values (Pa , Pb , Px , Pc . Pd ) of the pixel (Px) and the horizontally-neighboring 
pixels (Pa,Pb,Pc,Pd) for providing at an ouput (FOUT) said horizontally- 
filtered value (Out) of the pixel (Px). 
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"Post-processing method for reducing artifacts in block-coded digital 
images, and post-processing device for actuating such method." 



***** 



ABSTRACT 

5 A post-processing method for reducing artifacts in block-coded digital 

images, characterized by providing for: dividing an input image (In) into a 
plurality of image blocks (IB); for each image block (IB), estimating global 
features (GF) of said image block providing information on an average 
content of image edges along the horizontal and vertical directions of said 

10 image block (IB); for each pixel (Px) of an image block (IB) under 

examination, estimating local features (LF) for said pixel (Px) providing 
information on the content of image edges along the horizontal and vertical 
directions of an image area near said pixel (Px); modifying the value of said 
pixel (Px) according to both said global features (GF) of the image block 

15 (IB) to which said pixel (Px) belongs and said local features (LF) of the 

image area near said pixel (Px). (Figure 1). 
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